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ABSTRACT 

The oxygen and nitrogen abundance evolutions with redshift and galaxy stel- 
lar mass in emission-line galaxies from the Sloan Digital Sky Survey (SDSS) are 
investigated. This is the first such study for nitrogen abundances, and it pro- 
vides an additional constraint for the study of the chemical evolution of galaxies. 
We have devised a criterion to recognize and exclude from consideration active 
galactic nuclei (AGNs) and star-forming galaxies with large errors in the line flux 
measurements. To select star-forming galaxies with accurate line fluxes measure- 
ments, we require that, for each galaxy, the nitrogen abundances derived with 
various calibrations based on different emission lines agree. Using this selection 
criterion, subsamples of star-forming SDSS galaxies have been extracted from 
catalogs of the MPA/JHU group. We found that the galaxies of highest masses, 
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those with masses > 1O 1L2 M , have not been enriched in both oxygen and nitro- 
gen over the last ~ 3 Gyr: they have formed their stars in the so distant past that 
these have returned their nucleosynthesis products to the interstellar medium be- 
fore z=0.25. The galaxies in the mass range from ~ i0 1LO M Q to ~ 10 11,2 M Q do 
not show an appreciable enrichment in oxygen, but do show some enrichment 
in nitrogen: they also formed their stars before z=0.25 but later in comparison 
to the galaxies of highest masses; these stars have not returned nitrogen to the 
interstellar medium before 2=0.25 because they have not had enough time to 
evolve. This suggests that stars with lifetimes of 2-3 Gyr, in the 1.5-2 M Q mass 
range, contribute to the nitrogen production. Finally, galaxies with masses < 
10 11 M Q show enrichment in both oxygen and nitrogen during the last 3 Gyr: 
they have undergone appreciable star formation and have converted up to ~20% 
of their mass into stars over this period. Both oxygen and nitrogen enrichments 
increase with decreasing galaxy stellar mass in the mass range from ~lO n M to 
~10 10 M Q , then slightly decrease with further decrease of galaxy mass. 

Subject headings: galaxies: abundances - galaxies: evolution - galaxies: ISM: - 
Hll regions 



1. INTRODUCTION 



The study of various correlations between the global properties of galaxies is very im- 
portant for the understanding of their formation and evolution. Among th ese correlations, 



the on e between galaxy mass and metallicity is one of the most important. iLequeux et al. 



(119791 ) first showed that the oxygen abundance correlates well with total mass for a sam- 
ple of dwarf irregular and blue compact galaxies, in the sense that the higher the total 
mass, the higher the heavy element content. Since the galaxy mass is a poorly known 
parameter, the luminosity - metallicity relation h as of ten been considered instead of the 
mass - metallicity relation. iGarnett &: Shields I (119871 ) found that abundances of spiral 



nearby galaxies has been consic 


ered subsequently by many authors ( < 


Harnett & Shields 


1987 


Skillman, Kennicutt & Hodee 


1989: Vila-Costas & Ed 


munds II 1992 


Richer & McCall 


1995 


Pilvusdn 20011 


r, Zaritskv, Kennicutt & Huchra 1994 


: Pilvuein. Vflchez. & Contini 


2004 


Guseva et al. h 


2009, among others). 



In recent years, the number of available spectra of emission-line galaxies has increased 
dramatically due to the completion of several large spectral surveys. Measurements of emis- 
sion lines in those spectra have been carried out for abundance determinations and inves- 
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tigation of the luminosity - metallicity relation. Thus, iMelbourne &: Salzer I (120021 ) have 
considered the luminosity - metallicity relation for 51 9 galaxies in the KPNO International 
Spectro scopic Survey (K I SS). lLamareille et al. I ( 120041 ) for 6,387 galaxies in the 2dF Galaxy 



Survey, iTremonti et al. I (120041) for about 53,000 galaxies in the Sloan Digital Sky Survey 



(SDSS), and lAsari et al. I (120071 ) for 82,302 SDSS galaxies. 



It has been found, however, that, when mass can be determined, the mass-metallicity 
correlation is considerably tighter than the luminosity-metallicity correlation, suggesting 
that mass may be a more meaningful physical parameter than luminosity. In the last 



many investigators ( 


Li 


lv. Carollo. & Stockton 


2003; 


Savaglio et al. 


hood: 


Erb et al. 


2006; 


Cowie & Barker 


2008; 


Maiolino et al. 1 


2008; 


Lamareille et al. 


2009; 


Lara-Lopez et al. 


2009) 



In those investigations, oxygen abundances have been derived using various methods. The 
general conclusion from those studies is that the oxygen abundance change of star-forming 
galaxies over the last half of the age of the universe appears to be somewhat moderate, 
with A(log(0/H)) ~ 0.3 or lower. This change is comparable to or less than the scatter in 
the observed oxygen a bundances (s ee Fig. 10 of iLillv. Carollo. fc Stockton I (120031 ). Fig. 13 of 
Savaglio et al. I ( 120051 ) . and Fig.6 of lLara-Lopez et al. I (120091 )). 



Until now, no attention has been paid to the redshift evolution of nitrogen abundances 
in galaxies, despite the fact that they present several advantages for the study of the chem- 
ical evolution of galaxies. First, since at 12+log(0/H) > 8.3, secondary nitrogen becomes 
dominant and the nitrogen abundan ce increases at a faster rate than the oxygen abundance 
(IHenry. Edmunds fc Koppen 1120001 ). then the change in nitrogen abundances with redshift 
should show a larger amplitude in comparison to oxygen abundances and, as a consequence, 
should be easier to detect. Furt hermore, there is a time delay in the nitrogen production as 



comp ared to oxygen production (IMaeder Ill992l ; Ivan den Hoek fc Groenewegen 1119971 ; iPagel 



19971 ). This provides an additional constraint on the chemical evolution of galaxies. These 
reasons have led us to consider here not only the redshift evolution of oxygen abundances 
but also that of nitrogen abundances. To carry out such an investigation, it is necessary to 
derive accurate oxygen and nitrogen abundances. 

The emission line properties of photoionized H II regions are governed by its heavy ele- 
ment content and by the electron temperature distribution within the photoionized nebula. 
In turn, the latter is controlled by the ionizing star cluster spectral energy distribution and by 
the chemical composition of the H II region. The evolution of a giant extragalactic H II region 
associated with a star cluster is thus caused by a gradual change in time of the integrated 
stellar energy dist ribution due to stellar evolution. This has been the subject of numer- 



ous investigations ( jStasiriska Ill978l Il980i : iMcCall. Rybski fc Shields Ill985t iDopita fc Evans 
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19861 : iMov. Rocca-Volmerrange fc Fioc Il200ll ; IStasinska fc Izotov 1120031 : iDopita et al. 1120061 
among others). The general conclusion from those studies is that Hn regions ionized by 
star clusters form a well-defined fundamental sequence in different emission-line diagnostic 
diagrams. The existence of such a fundamental sequence forms the basis of various investi- 
gations of extragalactic H II regions. 



First, Baldwin. Phillips fc Terlevich I ( 119811 ) have suggested that the position of an ob- 
ject in some well-chosen emission-line diagrams can be used to separate H II regions ionized 
by star clusters from other types of emission-line objects. This idea has found general accep- 
tance and use. Thus, the [OIII]A5007/H/3 vs [NII]A6584/Ha diagram has been used widely 
to distinguish between star- forming galaxies and active galactic nuclei (AGNs). In partic- 
ular, the SPSS emission-line gal axies occupy a well-defined region shaped like the wings of 



a seagull (IStasinska et al. 1 120081 ) . The left wing consists of star-forming galaxies while the 



right wing is attributed to AGNs. However, t he exact location of t he dividing line between 



Hll regions and AGN s is still controversial (IKewley et al. I l2001t iKauffmann et al. 1 12003 



Stasihska et al. 1120061 ) . 



Second, iPagel et al. I (Il979l ) and lAUoin et al. I (119791 ) have suggested that the posi 



tions of Hll regions in some emission-line diagrams can be calibrated in terms of their 
oxygen abundances. This approach to abundance determination in Hll regions, usually re- 
ferred to as the "strong line method" has been widely adopted, especially in cases where 
the temperature-sensitive [O ill] A4363 line is not detected. Numerous relations have been 



ture estimates (' 


DoDita & Evans 


1986; 


Zaritskv. Kennicutt & Huchra 


1994; 


Pilvugin 


2000. 


2001c 


; IPilvugin & Thuan 


2005; 


Pettini & Pagel 


2004; 


Tremonti et al. 


2004; 


Liang et al. 


2006; 


Stasihska 


2006 


, among many others) 



It should be noted that the classic T e method, which relies on the electron temperature 
T e determined from the [O ill] A4363 line, is also based on the existence of the fundamental 
sequence of Hll regions. Indeed, when only one electron temperature, t3 or t2, is measured 
(t 3 and t 2 being the electron temperatures in the [O ill] and [O 11] zones, respectively) it is 
standard practice to use a t 2 - t 3 relation to estimate the other temperature. This relation is 
usually established on the basis of H II regions models which belong to the same fundamental 
sequence. 

Our study will be based on the SDSS data base of a million spectra. To study the redshift 
evolution of oxygen and nitrogen abundances, accurate oxygen and nitrogen abundance 
determinations are mandatory. The determination of accurate abundances in Hll regions 
from SDSS spectra poses two problems. First, line fluxes in SDSS spectra are measured 
by an automatic procedure. This inevitably introduces large flux errors for some objects. 
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We need to devise a way to recognize those objects and exclude them from consideration. 
Second, the SDSS galaxy spectra span a large range of redshifts. There is thus an aperture- 
redshift effect in SDSS spectra. Indeed, they are obtained with 3-arcsec-diameter fibers. At 
a redshift of z=0.05 the projected aperture diameter is ~ 3 kpc, while it is ~ 15 kpc at a 
redshift of ,2=0.25. This means that, at large redshifts, SDSS spectra are closer to global 
spectra of whole galaxies, i.e. to that of composite nebulae including multiple star clusters, 
rather than to spectra of individual Hll regions. One should then expect that some SDSS 
objects will not follow the fundamental H II region sequence. These need also to be identified 
and excluded from our sample. 

We propose here a method to recognize objects that suffer from one or both of these 
problems. It is based on the idea that if i) an object belongs to the fundamental Hll region 
sequence, and ii) its line fluxes are measured accurately, then the different relations between 
the line fluxes and the physical characteristic of H II regions, based on different emission lines, 
should yield similar physical characteristics (such as electron temperatures and abundances) 
of that object. 

The relations used for the determination of oxygen and nitrogen abundances and for 
the selection of star-forming galaxies with accurate line flux measurements are derived in 
Section 2. The selected subsamples of SDSS galaxies are described in Section 3. The relations 
between redshift, galaxy mass and metallicity are discussed in Section 4. Section 5 presents 
the conclusions. 

Throughout the paper, we will be using the following notations for the line fluxes: 
R 2 = [O n]A3727+A3729 = /[o//]A3727+A3729/-^H,a, 

N 2 = [Nll]A6548+A6584 = /[7V//]A6548+A6584//h / 3, 

R 3 = [Om]A4959+A5007 = / [O iii]A4959+A5007//h/3, 

The electron temperatures will be given in units of 10 4 K, and the stellar masses of galaxies 
in solar units. 



2. METALLICITY AND TEMPERATURE CALIBRATIONS 



We noted above that many relations have been proposed to convert various metallicity- 
sensitive emission-line combinations into metallicity or temperature estimates. The oxygen 
abundances derived with these different calibrations are not in good agreement, with differ- 
ences amounting up to 0.7 dex (IKewley fc Ellison II2008I ). In fact, there appears to exist as 
many different oxygen abundance scales as there are calibrations. Which metallicity scale 
is the correct one? There is strong evidence in favour of the metallicity scale defined by 
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the classic T P method, con sidered the most reliable one (see the discussion in iPilyugin 
20031 ; iBresolin et al. Il2009bl ). Indeed, the oxygen abundances derived with the T e method 
have been confirmed by high-precision model-independent determinations of the interstel- 
lar oxygen abundance in the solar vicinity, using high-resolution observations of the weak 
interstellar O IA1356 absorption line towards stars, and by recent determinations of stellar 
abundances. To put our derived abundances on the T e scale, we will use Hll regions with 
T e -based abundances as the basis of our calibrations. 

To establish the relations between different line fluxes and the physical characteristic of 
H II regi ons, we will be us i ng th e data for high-metallicity H II regions in nearby galaxies com- 
piled by IPilyugin et al. I (120091). We have added to this compilat i on the rec ent spectropho- 
tometric measurements of Saviane et al. I (b()08h : IBresolin et al~l (bo09afl ): Esteban et al 



(120091 ) . The oxygen and nitrogen abundances in those Hll regions were rederived in an 
uniform way within the framework of the standard Hll region model with a two- zone tem- 
perature distribution within the nebula. The equations f or th e electron temperatures and 
oxygen ionic abundances were taken from llzotov et al. I (120061 ). The relation between the 
electron temperature t2 within the + zone and t3, that within the ++ zone, is 



t 2 = 0.672*3 + 0.314. 



(1) 



( jCampbell. Terlevich fc Melnick 



1986 



This relation derived in IPilyugin et al. I (120091) is ve ry close to the commonly used one 



Garnett Ill9921 ). 



The abundance of the nitrogen ion N + is derived from the equation 

iV + N 2 0.726 
log — = log -± + 0.273 - — 0.02 log t 2 + 0.007£ 2 . 

The total nitrogen abundance is determined from 

.N.O.N 
log jj = log — + log — 



with the assumption 



7V+ 
0+ 



N 
O' 



(2) 



(3) 



(4) 



It is common practice, in constructing the calibration, to establish relations between the 
oxygen abundances and the strong-line fluxes. As a result, oxygen abundances are the best 
studied. In particular, the redshift evolution of oxygen abundances has been the subject of 
several recent investigations, as noted above. Here, we will adopt a different method. We 
will calibrate the positions of H II regions in three different diagnostic diagrams in terms of 
nitrogen abundances, to which less attention has been paid. The oxygen abundance will 
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then be estimated from the obtained nitrogen abundance and the N/O ratio using equations 
([2]) and ([3]). The independent determination of nitrogen abundances and of their redshift 
evolution will provide an independent check of the conclusions derived from the study of the 
redshift evolution of oxygen abundances. 

FigJTJ shows the nitrogen abundance as a function of log(N 2 ) (upper panel), log(N 2 /R 2 ) 
(middle panel), and log(N 2 /R 3 ) (lower panel). The filled circles are Hll regions with mea- 
sured electron temperatures ts^ (the auroral line [OIII]A4363 is available). The open circles 
are Hll regions with measured electron temperatures t2,jv (the auroral line [NII]A5755 is 
available). Those data have been fitted by the cubic expression 

12 + log(N/H) = a + axX + a 2 X 2 + a 3 X 3 (5) 

where X is successively X=log(N 2 ), X=log(N 2 /R 2 ), and X=log(N 2 /R 3 ). The derived curves 
are shown in FigJTjby solid lines. The constant coefficients in Eq.(jS]) for all three cases are 
given in Table [TJ 

FigfSJshows the electron temperature t 2 as a function of log(N 2 ) (upper panel), log(N 2 /R 2 ) 
(middle panel), and log(N 2 /R 3 ) (lower panel) for the same sample of Hll regions. Those 
data have also been fitted by the cubic expression 

t 2 = a + ai X + a 2 X 2 + a 3 X 3 (6) 

where X is successively X=log(N 2 ), X=log(N 2 /R 2 ), and X=log(N 2 /R 3 ). The derived curves 
are shown in Figj2]by solid lines. The constant coefficients in Eq.([6]) for all three cases are 
also given in Table [TJ 

The relations between nitrogen abundances and the abundance-sensitive indexes N 2 , 
N 2 /R 2 , and N 2 /R 3 will be used for the determination of nitrogen abundances in the SDSS 
galaxies. They will also serve to select star-forming galaxies with accurate line fluxes mea- 
surements. As for the relations between the electron temperature t 2 and the N 2 , N 2 /R 2 , 
and N 2 /R 3 indexes, they will be used for estimating the N/O ratios and deriving the O/H 
oxygen abundances from the N/H nitrogen abundances. 



SAMPLE SELECTION 



Line flux measurements in SDSS spectra have been carried out by several groups. We use 
here the data in several catalogs made available publicly by the MPA/JHU group Q. These 



The catalogs are available at 



http://www.mpa-garching.mpg.de/SDSS/ 
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catalogs give line flux measurements, redshifts and various other derived physical properties 
such as stellar masses for a lar ge sample of SPSS galaxies . The techniques used to construct 



the catalogues are described in iBrinchmann et al. I (12004J ) ; iTremonti et al. I ( 120041 ) and other 



publications of those authors. We have chosen to use these catalogues instead of the original 
SDSS spectral database because they contain generally more acc urate line flux measuremen ts 



[see a discussion of the errors in the line flux measurements by IBrinchmann et al. I (120041 ) ). 



In a first step, we extract from the MPA/JHU catalogs all emission-line objects with 
measured fluxes in the HB, Ha, [OII]AA3727,3729, [OIII]A4959, [OIII]A5007, [NII]A6548, 
[NII]A6584, [SII] A6717, and [SII] A6731 emission lines. The hydrogen, oxygen and nitrogen 
lines will serve to estimate oxygen and nitrogen abundances relative to hydrogen, and the 
ratio of the sulfur line intensities is an indicator of electron density. Since our calibrations 
are only valid in the low-density limit, we have included only those objects with a reasonable 
value of the [SII] ratio, i.e. those with 1.25 < F[sii]X67i7/F[Sii]X673i < 1-5. This results in 
a sample containing 118,544 objects which will be referred to hereafter as the total sample. 
The wavelength range of the SDSS spectra is 3800 - 9300 A so that for nearby galaxies with 
redshift z < 0.02, the [O n]A3727+A3729 emission line is out of that range. The absence of 
this line prevents the determination of the oxygen abundance, so all SDSS galaxies with z 
< 0.02 were also excluded. Thus, all galaxies in our total sample have redshifts greater than 
~ 0.023, i.e. they are more distant than ~ 100 Mpc. The redshift z and stellar mass Ms of 
each galaxy were also taken from the MPA/JHU catalogs. 

The emission-line fluxes are then corrected for interstellar reddening using the theoretical 
Ha to H B ratio and the analytical ap p roxim ation to the Whitford interstellar reddening 



law from llzotov. Thuan fc Lipovetsky I (119941 ) . In several cases, the derived value of the 



extinction c(H3) is negative and has been set to zero. 

For each galaxy, we have estimated three values of the nitrogen abundance and three 
values of the electron temperature t 2 , using the calibrations discussed above and the dered- 
dened N 2 , R 2 , R3 line intensities. Since measurements of the [NII]A6584 line are more reliable 
than those of the [NII]A6548 line, we have used N 2 = 1.33[NII]A6584 instead of the standard 
N 2 = [NII]A6548 + [NII]A6584. Then the mean value N/H of the nitrogen abundance for 
each galaxy is determined as 

, (AT/m log{N/H) N2 + \og{N/H) N2/R2 + \og{N/H) N2/R3 

log(iv / ti ) mean = . (7 J 

The mean value of the electron temperature t 2 is determined in a similar way. Using N/H 
and the mean t 2 , the N/O ratio and the oxygen abundance O/H are then derived using 
equations (J2]) and (J3J). 

We have suggested above that if i) the object belongs to the fundamental sequence of HII 
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regions, and ii) its line fluxes are measured accurately, then the calibrations based on different 
lines should result in similar abundances. We have calculated for each galaxy the deviations 
of individual values of the nitrogen abundance from the mean value Alog(N/H) n 2 ,N2/R2,n 2 /R3 
= log(N/R) N2tN2 / R2jN2/R3 - log(N/H) metm . The mean deviation Alog(N/H) mecm and the 
maximum deviation Alog(N/H) ma;r were also computed for each object. We have then used 
the value of Alog(N/H) maa; as a selection criterion to extract from the total sample three 
subsamples of star-forming galaxies with accurate line flux measurements, going from the 
most stringent requirement to the least stringent one : 

- Subsample A contains only objects with Alog(N/H) ma:r < 0.05, a total of 15,548 galaxies. 

- Subsample B contains only objects with Alog(N/H) maa , < 0.10, a total of 55,189 galaxies. 

- Subsample C contains only objects with Alog(N/H) max < 0.15, a total of 84,364 galaxies. 

Does selecting objects in such a way give us only 1) star-forming galaxies and 2) with 
accurate line flux measurements? We can test the reliability of our selection criterion con- 
cerning the first point by appealing to the [OIII]A5007/H/3 vs [NII]A6584/Ha diagram which 
is often used to distinguish between star-forming galaxies and AGNs. Figj3] shows such 
diagrams for subsamples A, B, C and for the total sample. The solid line, taken from 



Kauffmann et al. I (120031 ). shows the dividing line betweeen Hn regions ionized by star clus- 
ters and AGNs. Figj3] clearly shows that, while the total sample contains both star-forming 
galaxies and AGNs (lower right panel), our criterion does select out subsamples containing 
only star-forming galaxies. 

Concerning the second point, FigJH shows the classical R3 - R2 diagram for subsamples 
A, B, and C and the total sample. Each galaxy is plotted as a gray (light blue in the electronic 
version) open circle. For comparison, H II regions in nearby galaxies with ac curate line flux 



measurements, from the compilation of iPilyugin. Vflchez. fc Contini I (120041 ). are shown as 



black filled circles. Inspection of Fig 0] shows that the selected subsamples of galaxies occupy 
the same area in the R3 - R2 diagram as the Hn regions in nearby galaxies with accurate 
measurements, while the total sample covers a considerably larger area. Evidently, our 
selection criterion picks out objects which have line intensities that are in agreement with 
the well-measured line intensities of HII regions in nearby galaxies. 

We next investigate the O/H - N/O diagram (FigJ5J). As before, the SDSS galaxies 
from subsamples A, B, and C are shown as gray (light-blue in the electronic version) open 
circles, and th e H II regions in nearby galaxies with r ecent precise measurements of electron 



temperatures (jBresolin I l2007t iBresolin et al. Il2009al ) are shown as black filled circles. We 



see that the selected SDSS galaxies and the Hn regions in nearby galaxies with precise 
measurements lie in the same general region in the O/H - N/O diagram. However, the 
selected SDSS galaxies are shifted systematically toward lower N/O ratios. This slight 
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systematic shift (of ~ 0.2-0.3 dex in log N/O) can be understood in the following way. 
The selection criteria pick out galaxies with strong ongoing star formation, resulting in a 
higher star formation rate in the selected SDSS galaxies as compared to normal galaxies 
with the same mass. It has long been known that such a star formation burst would cause 
a temporary decrease of the N/O ratio in the galaxy (e.g. iPilyugin Ill992l . 119931 ). 

We now compare the properties of the galaxies in subsamples A, B, and C. The his- 
tograms of galaxy stellar masses, redshifts, oxygen and nitrogen abundances for the three 
subsamples are shown in FigJ6J The locations of the galaxies in the three subsamples in the 
redshift - galaxy stellar mass diagram are presented in FigJTl Examination of the distribu- 
tions of the parameters in the three subsamples reveals that, while they are somewhat similar 
for subsamples B and C, the distribution for subsample A differs significantly. In particular, 
there is a shift toward higher masses (by about 0.2 dex) of the galaxy mass distribution 
in subsample A relative to the other two subsamples. There is also a shift toward higher 
nitrogen abundances (by about 0.1 dex) of galaxies in subsample A relative to subsamples B 
and C (see Figj6]). These shifts seem to be caused by a too restrictive selection criterion. The 
selection condition for subsample A, Alog(N/H) max < 0.05, appears to be too constraining, 
eliminating too many galaxies and, as a consequence, causing a selection effect. In order not 
to bias our results, we will consider, in the remainder of the paper, subsample B as the basic 
subsample, while subsample C will be used as a control subsample. 



4. THE REDSHIFT EVOLUTION OF OXYGEN AND NITROGEN 

ABUNDANCES 

We investigate here the changes in oxygen and nitrogen abundances with galaxy stellar 
mass and redshift. Examination of Figj6]and FigJT] shows that such a study is justified only 
in the range of stellar galaxy masses from ~ 1O 9,5 M0 to ~ 10 1L5 M Q , where the number of 
galaxies is large enough to give good statistics. 



4.1. The Z-M5-O/H relation 



We first discuss the evolution of oxygen abundances with galaxy stellar masses. The 
upper panel of Fig JH] shows the oxygen abundances of galaxies in subsample B, with redshifts 
in the range 0.04 < z < 0.06, as a function of galaxy stellar mass. The lower panel of FigJH] 
shows the same diagram, but for galaxies in a higher redshift range, 0. 23 < z < 0.27. We note 
that our M5-O/H diagram is very similar to that of lErb et al. I ( 120061 ) (see their Fig.3). This 
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may not be surprising since those authors also used a N2 calibration to estimate the oxygen 
abundances of their SDSS galaxies. FigJH] shows that the oxygen abundance increases with 
increasing galaxy stellar mass up to a value M* s . For galaxies with M5 > M* s , the oxygen 
abundance becomes constant. Comparison of the upper and lower panels of FigJH] reveals 
that the value of depends on redshift, shifting to higher values at higher redshifts. 

To be more quantitative, we approximate the change in oxygen abundance with redshift 
and galaxy stellar mass by the following "redshift - galaxy stellar mass - oxygen abundance" 
z-Ms-O/H relation: 



12 + log(0/#) = a, + a 2 z + (a 3 + a 4 z) log(^), M s < M. 
12 + log(0/#) = ai + a 2 z, M s > M. 

log(M*/M ) = a 5 + a 6 z 



Fitting the data for galaxies in subsample B gives the following: 

12 + log(0/#) = 8.67 - 0.027 z + (0.41 -0.76 z)log(fff), M s < M* 

12 + log(0/#) = 8.67 - 0.027 z, M s > M* (9) 

log(M|/M ) = 9.60 + 5.65 z 

The derived z-M^-O/H relation for subsample B is shown in FigJH] by the solid line for 
2=0.05, and by the dashed line for z=0.25. 

To test the robustness of the derived Z-M5-O/H relation, we have also analyzed sub- 
sample C in a similar way. The data for galaxies in subsample C is shown in Figj9] Fitting 
the subsample C data gives the following: 

12 + log(0/#) = 8.67 - 0.020 z + (0.39 - 0.63 z) log(fff), M s < M* s 

12 + log(0/#) = 8.67 - 0.020 z, M s > M* s (10) 

log(M*/M ) = 9.54 + 5.90 z 

The Z-M5-O/H relation for subsample C is shown in Figj9] by the solid line for z=0.05, 
and by the dashed line for £=0.25. Comparison between FigJHJ and Figj9] shows that the 
obtained z-Mg-O/H relations for subsamples B and C are very similar, so that our results 
appear to be robust. The derived relations should approximate well the redshift evolution of 
oxygen abundances for z < 0.25, corresponding to lookback times of up to ~ 3 Gyr. Beyond 
z = 0.25, the number of galaxies in each subsample becomes very small (see FigJTJ) and 
the redshift evolution of abundances cannot be studied with enough statistics. To discuss 
abundance evolution, we will therefore restrict ourselves to objects with z < 0.25. At the 
lo w- redshift end, to reduce errors from aperture effects and following the recommendations 



of iKewley. Jansen &: Geller I (120051 ). we will limit ourselves to objects with redshifts z > 0.04. 



Objects with 0.04 < z < 0.06 will be considered as representative of the present-day epoch. 
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FigJTU] shows the distribution of the deviations of the oxygen abundances from the 
derived z-M^-O/H relation for both subsamples B and C. It is seen that the scatter of oxygen 
abundances about the Z-M5-O/H relation is slightly larger for galaxies in subsample C than 
in subsample B. This effect can also be seen by comparing FigJH]to FigJHJ Nevertheless, it 
can be said that, for both subsamples of galaxies, the oxygen abundances follow reasonably 
well the derived z-M^-O/H relation. The mean deviation is small, being ~ 0.06 dex. We 
thus conclude that the derived relations can be used, in principle, to obtain a rough estimate 
of the metallicities of SDSS galaxies. 



4.2. The metallicity plateau and the oxygen enrichment as a function of 

galaxy stellar mass 

One of the most remarkable features of the M^-O/H relation at the current epoch is the 
metallicity plateau at high galaxy stellar masses. For the sake of defmiteness, we will discuss 
the results for subsample B. The oxygen abundance increases with galaxy stellar mass until 
Ms = 1O 9,9 M , but then remains approximately constant, equal to 12+log(0/H)=8.67, in 
galaxies with M s > 10 9 - 9 M Q (see the upper panel of FigJS]). What is the physical meaning of 
such a plateau? The observed oxygen abundance in a galaxy is defined by its astration level 
or gas mass fraction /i, and by the mass exchange between a galaxy and its environment 
(ISearle fc Sargent Ill972l ; IPagel 1119971 ). It is believed that ga lactic winds do not play a signif- 



icant role in the chemi cal evolution of large spiral galaxies ( jGarnett 1 12002| ; iTremonti et al. 



2004: Dalcanton 



2007 ) and that the rate of gas infall onto the disk de creases exponentiall y 



with time ( Matteucci fc Francois 1989 ; Pilvugin fc Edmunds 1996al lbl; Icalura et al. 2009 ). 
The present-day location of a system in the /i-O/H diagram is then governed by its evolu- 
tion in the recent past, a nd is only weakly dependent of its evolution on long time-scales 
(jPilyugin fc Ferrini Ill998l ) . Therefore the observed oxygen abundance in a large spiral galaxy 
is mainly defined by its gas mass fraction /1. The presence of a metallicity plateau in the 
Ms-O/H diagram at all redshifts implies then that the high-mass SDSS galaxies have very 
similar gas mass fractions. 

The luminosity- central metallicity relation for near by spiral galaxies also shows a plateau 
at high luminosities (jPilyugin. Thuan fc Vflchez 1120071 ) . This plateau was interpreted as ev- 
idence that the gas in the centers of the most metal-rich galaxies has been almost completely 
converted into stars and that the oxygen abundance in the centers of the most luminous 
metal-rich galaxies has reached its maximum attainable value of 12+log(0/H) ~ 8.87. The 
plateau for oxygen abundances in the SDSS galaxies is at the level of 12+log(0/H)=8.67, 
i.e. 0.2 dex smaller. The simple model of chemical evolution of galaxies predicts that a 



13 



decrease of /i by 0.1 results in an increase of oxygen abunda nce by ~0.13 dex, in the range 



of \x from ~0.50 to ~ 0.05 (jPilyugin. Thuan &: Vflchez 1120071 ). Then, the difference between 



the maximum attainable value of the oxygen abundance and the mean oxygen abundance 
in high-mass SDSS galaxies, A(log(0/H)) ~ 0.2, corresponds to a difference in gas mass 
fraction A/i ~ 0.15. Since the maximum attainable oxygen value corresponds to complete 
astration, i.e. \i = 0, the mean /i in high-mass SDSS galaxies at the present epoch should 
be ~ 15%, with a probable range of /i from ~ 5% to ~ 25%. 

Why do our SDSS subsamples not contain galaxies with abundances as high as the 
maximum attainable value? The reason has to do with the gas content of the SDSS galaxies 
in our subsamples. As said before, the maximum attainable value of the oxygen abundance 
corresponds to the limiting case where the gas has been completely converted into stars. 
The galaxies in our SDSS subsamples are all characterized by strong emission lines in their 
spectra, meaning that they are undergoing strong starbursts. This requires in turn that they 
contain an appreciable amount of gas. In other words, because our galaxies are gas-rich, 
they have not reached the maximum attainable value of the oxygen abundance which re- 
quires complete gas exhaustion. If galaxies were not gas-rich, they would have weak emission 
lines, unlikely to be measured accurately. There exists thus a lower limit on the gas mass 
fraction for a galaxy to have accurately measured lines and to be included in our subsam- 
ples. In that sense, the observed plateau is likely affected by the selection criteria. We 
note however that, even in the most evolved nearby spiral galaxies, the Hll region oxygen 
abundances are generally lower than the maximum attainable value. There is furthermore 
an aperture effect that slightly lowers the metallicities observed for SDSS galaxies. The 
maximum attainable value oxygen abundance observed in the central part of some nearby 
galaxies is usually derived by linearly fitting the variations of H II region abundances with 
galactocentric distance, and extrapolating to R=0. In the distant SDSS galaxies, because 
one fiber includes many H II regions that show decreasing metallicities towards larger galac- 
tocentric distances, the oxygen abundances are slightly diluted when averaged over a large 
region. 

In the lower panel of FigJHl we compare the M5-O/H relation for local (z ~ 0.05) galaxies 
(solid line) with that for distant (z « 0.25) ones (dashed line). Three features are to be noted. 
First, it can be seen that, for the galaxies of highest stellar masses, those with masses > 
10 11 M Q , the plateau value of the oxygen abundance does not change in the redshift interval 
from 0.05 to 0.25. This implies that the galaxies of highest mass in our SDSS subsamples 
have reached their highest astration level some 3 Gyr ago, and have been somewhat "lazy" 
in their evolution afterwards. Second, for lower mass galaxies with stellar masses < 10 11 M Q , 
it is seen that the value of the oxygen enrichment during the last 3 Gyr increases with 
decreasing galaxy mass, in the mass interval from 10 n M Q to 10 9 ' 9 M Q , as shown by the 



-14- 



widening gap between the solid and dashed curves toward lower masses. At Ms=lO 9 ' 9 M , 
there is a difference Alog(0/H) ~ 0.25 between a local galaxy and one at redshift 0.25. For 
galaxies with M s < 1O 9 - 9 M , the oxygen enrichment during the last 3 Gyr slightly decreases 
with decreasing galaxy mass, the slope of the solid line being slightly steeper than that of 
the dashed line. Third, the value of M* s , the mass where the oxygen abundance becomes 
constant with galaxy stellar mass, is redshift-dependent, becoming higher at larger redshifts. 

In summary, analysis of the z-M^-O/H relation has led to the following main conclu- 
sions. 

- The galaxies of highest masses, those with Ms > 1O 11 M , have reached their highest astra- 
tion level in the past and have not had an appreciable oxygen abundance enrichment during 
the last ~ 3 Gyr. 

- The mean value of the oxygen enrichment during the last 3 Gyr in galaxies with stellar 
masses in the range from 1O 1O M to 1O U M is A(log(0/H)) ~ 0.11, with A(log(0/H)) ~ 
0.23 at 1O 1O M and A(log(0/H)) =0 at 1O U M . 

The above picture will now be put to the test through the study of the redshift evolution 
of nitrogen abundances, which we consider next. 



4.3. The z-Ms-N/H relation 

We proceed in the same way as in our analysis of the oxygen abundances. The upper 
panel of Fig {11] shows the nitrogen abundances of the local galaxies in subsample B, with 
redshifts 0.04 < z < 0.06, as a function of galaxy stellar mass. The lower panel of FigfTTl 
shows the same diagram, but for more distant galaxies, with redshifts 0.23 < z < 0.27. FigfTTl 
shows that the general behavior of nitrogen abundances with redshift and galaxy stellar mass 
is similar to that of oxygen abundances. The nitrogen abundance increases with increasing 
galaxy stellar mass up to a value M* s of the stellar mass. Then, for galaxies with M5 > M* s , 
the nitrogen abundance remains approximatively constant, reaching a plateau. The value of 

is redshift-dependent, becoming higher at larger redshifts. 

The change of nitrogen abundance with redshift and galaxy stellar mass can be approx- 
imated by a Z-M5-N/H relation, similar to the one for oxygen abundances. Fitting the data 
for galaxies in subsample B results in the following relation: 

12 + log(N/H) = 7.89 -0.110 z +(0.99 - 1.40 z) log(§J), M s < M s 

12 + \og(N/H) = 7.89- 0.110 z, M s > M* s (11) 

\og(M* s /M & ) = 9.97 + 4.92^ 

The derived relation is shown in FigflTJby a solid line for z=0.05, and by a dashed line for 
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2=0.25. 

Again, to test the robustness of the obtained z-Mg-N/H relation, we have examined 
subsample C in a similar way. The data for galaxies in subsample C are plotted in FigJTSZl 
Fitting those data gives: 

12 + log(N/H) = 7.88 -0.027 z + (0.85 - 1.14 z) log(^), M s < M* 

12 + \og(N/H) = 7.88 - 0.027 z, M s > M* (12) 

log(M|/M ) = 10.04 + 5.03 z 

The derived z-Mg-N/H relation is shown in Fig{12] by a solid line for 2=0.05, and by a 
dashed line for 2=0.25. Comparison of Figures [TT1 and [121 shows that they are very similar 
and that the derived z-Mg-N/H relation is robust. 

While the general behavior of nitrogen abundances with redshift and galaxy stellar mass 
is quite similar to that of oxygen abundances (compare FigsJH]and[TT]), there are significant 
differences, due to the different production mechanisms of these two elements. The de- 
pendence of the nitrogen abundance on galaxy stellar mass for Mg < M* s is considerably 
steeper than that for oxygen abundances. This is caused by the fact that at oxygen abun- 
dances higher than about 12+log(0/H)=8.3, the metallicity-dependent nitrogen production 
by intermediate-mass stars starts to dominate. Another remarkable difference concerns M* s , 
the mass which marks the transition from the linear regime to the plateau regime: at all 
redshifts, M s is shifted towards higher values in the Mg-N/H diagram as compared to in 
the Mg-O/H diagram. Thus, logM 5 is equal to 10.2 and 11.2 at z = 0.05 and z = 0.25 
respectively in the Mg-N/H diagram, as compared to 9.9 and 11.0 in the Mg-O/H diagram. 

We now attempt to understand this M* s shift. The time delay in nitrogen production 
relative to oxygen production plays an important role. Examination of FigsIEl and [TT] shows 
that galaxies with masses > 10 1L2 M Q are in the plateau regime in both the Mg-O/H and 
in Mg-N/H diagrams, at both 2=0.05 and 2=0. 25. This means that those galaxies have not 
undergone appreciable enrichment in both oxygen and nitrogen during the last ~ 3 Gyr. 
Evidently, there has not been appreciable star formation in those galaxies over the redshift 
range from 2=0. 25 to ^=0.05. Significant star formation in those galaxies has occurred so 
long ago that stars have returned their nucleosynthesis products to the interstellar medium 
before the epoch corresponding to 2=0.25. 

Galaxies with masses between ~ 1O 1LO M and ~ 1O 1L2 M do not show an appreciable 
enrichment in oxygen abundance from 2=0.25 to 2=0.05 but do show some enrichment in 
nitrogen over this period. This suggests that there has not been appreciable star formation 
in those galaxies over the period from 2=0.25 to 2=0.05. However, they do contain stars 
that were formed before 2=0.25, but later in comparison to the galaxies of highest masses. 
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The massive oxygen-producing stars die after a few million years, releasing oxygen in the in- 
terstellar medium. By contrast, the nitrogen-producing intermediate-mass stars have longer 
lifetimes, and so they have not returned nitrogen to the interstellar medium before z=0.25 
because they have not had enough time to evolve. This also suggests that stars that make 
a contribution to the nitrogen production have lifetimes of a few Gyr. 

The galaxies with masses < 1O 1LO M show enrichment in both oxygen and nitrogen 
abundances after the period corresponding to z=0. 25. This means that appreciable star 
formation has taken place in those galaxies during the last ~ 3 Gyr. The nitrogen production 
increases with decreasing galaxy mass from 1O 1L2 M to ~ 10 ia2 M Q where it reaches a value 
Alog(N/H) ~ 0.65, then it decreases with further decrease of galaxy mass. 

The lower panel of FigJTUl shows the deviations of nitrogen abundances from the derived 
z-Mg-N/H relation for galaxies from both subsamples B (solid line) and C (dashed line). 
Examination of the upper and lower panels of FigJTU] shows that the lower histogram is 
broader than the upper one, i.e. that the nitrogen abundances show a larger scatter around 
the z-Mg-N/H relation as compared to oxygen abundances around the z-Mg-O/H relation. 
The mean deviation for nitrogen abundances is A(log(N/H)) ~ 0.15 against A(log(OH)) ~ 
0.06 for oxygen abundances. However, nitrogen also spans a significantly larger abundance 
range, from 12+log(N/H) ~ 6.5 to 12+log(N/H) ~ 8.0, as compared to oxygen which goes 
only from 12+log(0/H) ~ 8.25 to 12+log(0/H) ~ 8.75. Taking into account the differences 
in range, then the relative scatters around the z-Mg-N/H and z-Mg-O/H relations are 
comparable. It should be emphasized that the larger scatter of nitrogen abundances around 
the z-Mg-N/H relation in comparison to the scatter of oxygen abundances around the z-Mg- 
O/H relation cannot be attributed to larger errors in nitrogen abundance determinations. 
Indeed, the O/H values are derived from the N/H and N/O values, and any error in the 
nitrogen abundance determination is propagated into the error in the oxygen abundance 
determination. Then, if the large scatter of nitrogen abundances around the z-Mg-N/H is 
caused by large errors in the nitrogen abundance determinations, then the oxygen abundances 
would show a similar or larger scatter around the z-Mg-O/H relation. Just the opposite is 
observed. 

The larger scatter of nitrogen abundances as compared to that of oxygen abundances 
can be understood in the following way. The oxygen abundance in a galaxy is mainly 
defined by its astration level. In contrast, the nitrogen abundance is defined not only by 
the astration level but depends also on the star formation history of the galaxy. As noted 
abov e, a star formation bu rst results in a temporary decrease of the N/O ratio in a galaxy 
(e.g. iPilyugin Ill992l Il993l ). As a result, galaxies with a given oxygen abundance can have 



different N/O ratios and nitrogen abundances, depending on their star formation histories. 
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We now check that expectation with our data. We use subsample C to have the most 
statistics. Fig{13] shows the oxygen and nitrogen abundances as a function of redshift for 
galaxies with masses in the lower range lO iao M -lO ia3 M (upper left panel) and in the 
upper range 1O 11,2 M -1O 1L5 M (upper right panel). The N/O abundance ratios of those 
galaxies as a function of redshift are shown in the lower panels. It can be seen that the 
galaxies in the two different mass ranges show different behaviors with redshift, reflecting 
the difference in the level of star formation activity in them. That in the lower mass range 
is high (those galaxies show an appreciable increase of oxygen abundance with decreasing 
redshift) while that in galaxies in the upper mass range is very low (those galaxies do not 
show a significant oxygen enrichment with decreasing redshift). FigflBl shows that, in spite 
of the different behaviors of the oxygen abundances with redshift in the two galaxy mass 
ranges, the scatter in oxygen abundances at a given redshift is comparable in the two cases, 
the dispersion of the data points in the low mass range being only slightly larger than that 
in the high mass range. By contrast, the scatter in the N/O ratios and in the nitrogen 
abundances at a given redshift is considerably higher for galaxies with high star formation 
activity (those in the low mass range, left upper and lower panels in Fig JT3|) than for galaxies 
with low star formation activity (those in the high mass range, right upper and lower panels 
in Fig JTBl . At z=0.05, oxygen abundances in galaxies with masses lO lo ' o M0-lO ia3 M are 
similar to those in galaxies with masses 1O 11 ' 2 M0-1O 11 ' 5 M . The maximum value of the 
N/O ratio in the two galaxy mass ranges is also similar. However, many galaxies in the 
low mass range possess low N/O ratios (FigJT3| lower left panel). The fact that the N/O 
ratios of many galaxies with high star formation activity are significantly lower than those in 
galaxies with low star formation activity does confirm our expectation that the larger scatter 
of nitrogen abundances around the z-Mg-N/H in comparison to that of oxygen abundances 
around the z-M^-O/H is caused by the temporary decrease of the N/O ratio, due to star 
formation bursts of different amplitudes and/or ages, i.e. to different star formation histories 
in galaxies. 

Thus, the consideration of the nitrogen abundance evolution with redshift and galaxy 
stellar mass has confirmed the general picture obtained from the oxygen abundance evolution 
analysis. Examination of both the z-Mg-N/H and the z-Mg-O/H relations has led to the 
following conclusions: 

- The galaxies of highest masses, those with masses > 1O 1L2 M , have reached their high 
astration level more than 3 Gyr ago, so that stars in those galaxies have returned their 
nucleosynthesis products to the interstellar medium before £=0.25. 

- The galaxies with masses in the range from ~ 1O 11O M to ~ 1O 1L2 M also form their stars 
before £=0.25, but later in comparison to the galaxies of highest masses. The intermediate- 
mass stars in those galaxies have not returned nitrogen to the interstellar medium before 
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z=0.25 because they have not had enough time to evolve. 

- Significant star formation has occurred in galaxies with masses lower than ~ 1O 11 M during 
the last 3 Gyr. Those galaxies have converted up to 20% of their total mass into stars over 
this period. 

- Stars with lifetimes of a few Gyr contribute to the nitrogen production. 



4.4. Influence of the selection criteria on the redshift — galaxy mass — 

metallicity relations 

We wish to examine here how our selection criteria affect the results obtained. We 
compare the results obtained for our SDSS subsamples with those obtained for a sample of 
star-forming SDSS g alaxies selected with commonly used criteria. We use the diagnostic 



diagram proposed by iBaldwin. Phillips fc Terlevich I ( 119811 ) where the excitation properties 



of Hn regions are studied by plotting the low-excitation [N ll]A6584/Ha line ratio against 
the high-excitation [O lll]A5007/H/3 line ratio. The diagram can be used to separate different 
types of emission-line objects according to their main excitation mechanism. We have thus 
excluded from the total SDSS sample all objects above the solid line given by the equation 

log([0 ,„]A5007/H 3 ) = hi[|Nl , 1 ^) Ha) . M5 + 1.3 (13) 



which separates objects with H II spectra from those containing an AGN (jKauffmann et al 



20031 ). We have extracted from the total sample a subsample of 95,046 star-forming galaxies, 



referred to hereafter as the "all star-forming galaxies" or SFG subsample. 

We compare the z-M^-O/H and z-Mg-N/H diagrams for the SFG and C subsamples. 
The upper panel of FigJTJ] shows the oxygen abundances of galaxies in the SFG subsample, 
with redshifts in the range 0.04 < z < 0.06, as a function of galaxy stellar mass. The lower 
panel of FigJT4l shows the same diagram, but for galaxies in a higher redshift range, 0.23 < z 
< 0.27. The derived z-Mg-O/H relations for subsample C, superimposed on the data points 
of the SFG subsample, are shown in FigJT4]by the solid line for z=0.05, and by the dashed 
line for z=0.25. Inspection of FigJT41 shows that the galaxies from the SFG subsample follow 
well the z-M^-O/H relation derived for subsample C. Comparison of FigJT4]with Figj9]shows 
however that the scatter of the points about the z-M^-O/H relation is larger for the SFG 
subsample than for subsample C. This is evidence that our selection criteria effectively weeds 
out galaxies with less reliable measurements. 

FigJT5l shows the Z-M5-N/H diagram. As before, the z-M^-N/H relations derived for 
subsample C are superimposed on the datapoints of the SFG subsample. Again, the galax- 
ies from the SFG subsample follow well the z-Mg-N/H relation derived for subsample C. 
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Thus, the results obtained with the commonly used method of selecting star-forming galax- 
ies do not differ appreciably from those obtained with our selection method. We emphasize 
that our approach presents however two important advantages. First, it is not necessary 
to know a priori the precise location of the dividing line between Hll regions and AGNs, a 



subject which is, as noted b efore, still controversial ( jKewley et al. Il200ll ; iKauffmann et al. 



20031 ; IStasihska et al. II2006I ). Second, our approach allows to reject all unreliable measure- 



ments. The counterpart of these advantages is that our approach requires a good sample of 
calibration datapoints - a reasonably large sample of H II regions with precise spectrophoto- 
metry measurements, including weak auroral lines - to establish reliable relations between 
the strong line fluxes and the physical characteristics of H II regions. Those relations should 
however be established in any case, as they are necessary for abundance determinations in 
star-forming galaxies, independently of the particular method used to select them. 



4.5. Comparison to previous work on the z— M5— O/H relation 

The evolution of the mass-metallicity relation of galaxies with redshift has been con- 
sidered previously by several groups, as described in the introduction. But, as discussed 
before, each group uses a different calibration to derive abundances which show as a result 
sometimes large discrepancies. So it is difficult to directly compare, or put on the same 
scale, the abundances derived by other groups and our own. Thus, we will not attempt 
such a comparison. Rather, we will limit ourselves to comparing the evolution in oxygen 
abundances with redshift which seems to be less sensitive to the adopted calibration. As for 
the nitrogen evolution with redshift, it has not been considered before. 



We first summarize the results obtained by previous investigators. iLilly. Carollo. fc Stockton 



( 120031 ) have estimated the oxygen abundance in a sample of 66 Canada-France Redshift Sur- 
vey galaxies in the redshift range 0.47 < z < 0.92, using the flux ratios of bright oxygen 
emission lines. They concluded that, at half the present age of the universe, the overall oxy- 
gen abundance of the galaxies in their sample, with luminosities ranging from M#=-20 to 
M#=-22 (or log(L#/L£o) ~ 10.2 - 11.0), is only slightly lower than the oxygen abundance 
in similar luminous galaxies today. They found a variation A(log(0/H)) = 0.08±0.06. 



Savaglio et al. I (120051 ) have investigated the mass-metallicity relation using galaxies 
from the Gemini Deep Deep Survey and the Canada-France Redshift Survey in the redshift 
range 0.4 < z < 1.0. Their galaxies with Ms > 1O 1O M have oxygen abundances close to 
those in local galaxies of comparable mass, while the oxygen abundances in galaxies with 
Ms < 1O 1O M are lower on average (with a large scatter) than those in galaxies of similar 
masses at the present epoch (see their Fig. 13). 
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Cowie h Barger I (120081 ) have studied the oxygen abundance evolution from z = 0.9 to z 
= 0.05 using a large sample of galaxies in the Great Observatories Origins Deep Survey-North 
(GOODS-N). They have found an evolution of the metallicity-mass relation corresponding 
to a decrease of 0.21±0.03 dex between the value at z — 0.77 and the local value in the 10 10 - 
1O 11 M range. They also found that star formation in the most massive galaxies (>1O U M ) 
ceases at 2<1.5 because of gas starvation. 



Lamareille et al. I (120091 ) have derived the mass-metallicity relation of star-forming galax- 
ies up to z ~ 0.9 using data from the VIMOS VLT Deep Survey. They found that the galaxies 
of 10 10 ' 2 solar masses show a larger oxygen enrichment (A(log(0/H)) ~ 0.28) from z ~ 0.77 
to z—0) than the galaxies of 10 9 ' 4 solar masses (A(log(0/H)) ~ 0.18). 



Lara-Lopez et al. I (120091 ) have studied the oxygen abundance of relatively massive 
(log(Ms/M Q ) > 10.5) star-forming galaxies from SDSS/DR5 at different redshift intervals 
from 0.4 to 0.04. They found an oxygen enrichment A(log(0/H)) ~ 0.1 from redshift 0.4 to 
0. 



Asari et al. I ( 120091 ) have derived the mass-metallicity relation at different lookback 
times for SDSS galaxies using the stellar metallicities estimated with their spectral synthesis 
code. They have found that the more massive galaxies show very little evolution since a 
lookback time of 9 Gyr. 

Examination of all these studies shows good qualitative agreement between them, but 
with a rather large scatter in the estimated values of the oxygen enrichment. The results 
of these investigations can be summarized as followed: 1) the most massive galaxies (those 
with masses > 1O 11 M ) do not show an appreciable enrichment in oxygen from z ~ 0.7 to z 
= 0; 2) in the 10 10 -10 11 M Q mass range, an increase of the oxygen abundance A(log(0/H)) 
~ 0.08-0.28 is observed in the redshift range from z ~ 0.7 to z ~ 0. Our results are also 
in good qualitative agreement with those previous results. We also see no change in oxygen 
abundance with redshift for galaxies with masses greater than 10 11 M Q , in the redshift range 
from ~0.25 to 0. We estimate the mean increase of the oxygen abundance with redshift in 
the 10 10 -10 n M o galaxy stellar mass range to be A(log(0/H)) ~ 0.11, with A(log(0/H)) ~ 
0.23 at 1O 1O M and A(log(0/H)) =0 at 10 11 M Q . This is in agreement with the upper range 
of values found in previous works, especially when we take into account the fact that we 
have considered a smaller redshift interval than previous investigators. 
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4.6. Nitrogen enrichment by intermediate mass stars: a discrepancy with 

current stellar evolution theory 



Edmunds &: Pagel I (119781 ) have suggested that observations of the N/0 abundance ratio 
in galaxies can be understood if N is manufactured in stars of 1-2.5 M . The N/0 ratio of a 
galaxy then becomes an indicator of the time that has elapsed si nce the bulk of star formation 



occurr ed, or in other words of the nominal "age" of the galaxy. iPilyugin. Thuan fc Vflchez 



(120031 ) have found that the N/O ratios in Hll regions of galaxies of early morphological types 
are systematically higher than those in Hll regions with the same O/H value in galaxies 
of late morphological types. Moreover, it is known that the star formation histories of 
galaxies of different morphological types differ, in the sense that the spiral galaxies of early 
morphological types hav e a significantly larger fraction of old stars than the galaxies of 



late morphological type (ISandage II1986I ). These two facts lead to the conclusion that the 
contribution to the nitrogen production by low- and intermediate-mass stars, with a long 
time delay of a few Gyr between their moment of birth and the time of release of their 
nucleosynthetic products in the interstellar medium, can be significant. 

One of the main conclusions of this study is that there has been a significant evolution 
in nitrogen abundances in star-forming galaxies over the last 2-3 Gyr. This appears not to 
agree with the predictions of current theoretical models of intermediate mass stars. Indeed, 
a time duration of 2-3 Gyr corresponds to the lifetime of stars in the 1.5-2 M mass 
range. While stellar evolution models of intermediate mass stars have long predicted that 
they contribute significantly to nitrogen enrichment, it is stars of 3 - 8 M Q that are supposed 
to do the job, not stars in the 1.5-2 M Q mass range. Intriguingly, there are other types 
of observations that also suggest that stars in the 1.5 - 2 M Q mass range contribute to 



nitrogen enrichment. iRicher fc McCall I (120081 ) have found that planetary nebulae in nearby 



galaxies (with or without ongoing star formation) often show large nitrogen enrichments. 
They have argued that these planetary nebulae are the descendants of relatively low mass 
progenitors, of approximately 1.5M Q or less, i.e. that low- and intermediate-mass stars are a 
more important source of nitrogen than has been hitherto considered. Thus, our conclusion 
that stars with lifetimes of a few Gyr contribute to the nitrogen production, derived from 
the consideration of the nitrogen abundance evolution with redshift in galaxies, is in line 
with the conclusions of other types of investigations. While stellar evolution theory does 
not yet predict nitrogen production in stars with masses of ~ 1.5 - 2 M , several lines of 
observations now do so. 
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5. CONCLUSIONS 

The redshift evolution of oxygen and nitrogen abundances in emission-line SDSS galaxies 
has been studied. 

We have paid particular attention to the construction of our galaxy sample, using the 
MPA/JHU catalogs of line flux measurements and other derived physical properties for 
SDSS galaxies. We have devised a way to recognize and exclude from consideration not only 
AGNs, but also star-forming galaxies with large errors in their line flux measurements. We 
have found that the requirement that nitrogen abundances, derived with different calibration 
relations based on different emission lines, agree, can be used as a reliable criterion to select 
star-forming galaxies with accurate line fluxes measurements. 

We have derived relations between nitrogen abundances and the abundance-sensitive 
N 2 , N 2 /R 2 , and N2/R3 indexes. Those relations have been used to determine nitrogen 
abundances in the SDSS galaxies. The small dispersion among the various derived nitrogen 
abundances for a given galaxy is used as a criterion to select star-forming galaxies with 
accurate line fluxes measurements. The relations between the electron temperature t 2 and 
the N 2 , N2/R2, and N2/R3 indexes have also been established. Those calibrations have been 
used to estimate the N/O ratio and derive the oxygen abundances O/H from the nitrogen 
abundances N/H. 

Subsamples of star-forming SDSS galaxies have been extracted from the MPA/JHU cat- 
alogs, using the small nitrogen abundance dispersion criterion described above. The nitrogen 
and oxygen abundances are estimated for these galaxies. The evolution of the oxygen and 
nitrogen abundances with redshift and galaxy stellar mass of galaxy are investigated, that 
of nitrogen abundances for the first time. We have obtained the following main results. 

1) The galaxies of highest masses (those more massive than ~ 1O 1L2 M ) do not show an 
appreciable enrichment in both oxygen and nitrogen from £=0.25 to ,2=0.05. Those galaxies 
have reached their high astration level in such a distant past that their stars have returned 
their nucleosynthesis products to the interstellar medium before £=0.25. 

2) The galaxies in the mass range from ~ 1O 1LO M0 to ~ 1O 1L2 M do not show an oxygen 
enrichment, but do show some enrichment in nitrogen. Those galaxies also formed stars 
before z=0.25, but at a later epoch in comparison to the galaxies of highest masses. Their 
stars have not returned nitrogen to the interstellar medium before £=0.25 because they have 
not had enough time to evolve. 

3) The galaxies with masses lower than ~ 10 11 M Q show enrichment in both oxygen and 
nitrogen abundances over the redshift period from £=0.25 to £=0.05, i.e. during the last 
3 Gyr. The oxygen enrichment increases with decreasing galaxy mass, from lO n M to 
M s =10 9 - 9 M o . It reaches a value Alog(0/H) ~ 0.25 at M 5 =10 9 - 9 M o and slightly decreases 
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with further decrease of galaxy mass. The nitrogen enrichment increases with decreasing 
galaxy mass, from ~ lO n M0 to ~ lO ia2 M . It reaches a value Alog(N/H) ~ 0.65 at ~ 
1O 1O ' 2 M and slightly decreases with further decrease of galaxy mass. Significant star forma- 
tion has occurred in those galaxies during the last 3 Gyr. They have converted up to 20% 
of their total mass to stars over this period. 

4) Stars with lifetimes of 2-3 Gyr, i.e. in the 1.5 - 2 M mass range, contribute to the 
nitrogen production. This is not in agreement with current stellar evolutionary models of 
intermediate mass stars which predict that stars in the 3-8 M mass range do the job, not 
stars in the 1.5 - 2 M mass range. 

5) The general picture of the oxygen abundance evolution with redshift and galaxy stellar 
mass obtained here and in previous work is confirmed and strengthened by consideration of 
the nitrogen abundance evolution. 
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Table 1. Values of coefficients in Eq.(jS]) for 12+log(N/H), and in Eq.(j6]) for t 2 



Y value 


X value 


a 




&2 


a 3 


12+log(N/H) 


log(N 2 ) 


7.649 


1.454 


0.257 


-0.148 


12+log(N/H) 


log(N 2 /R 2 ) 


7.918 


0.877 


-0.058 


0.038 


12+log(N/H) 


log(N 2 /R 3 ) 


7.526 


0.521 


0.062 


0.014 


t 2 


log(N 2 ) 


0.778 


-0.354 


0.001 


0.044 


h 


log(N 2 /R 2 ) 


0.702 


-0.248 


0.042 


0.019 


t 2 


log(N 2 /R 3 ) 


0.802 


-0.172 


0.003 


0.011 
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Fig. 1. — Nitrogen abundance as a function of different abundance-sensitive indexes for Hll 
regions in nearby galaxies with measured electron temperatures t^o or t2,jv- The filled circles 
show H II regions with measured t^o temperatures. The open circles show H II regions with 
measured t 2 ,jv temperatures. The solid lines are the adopted cubic fits. 
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Fig. 2. — Electron temperature t2 as a function of different abundance-sensitive indexes for 
Hn regions in nearby galaxies with measured electron temperatures t^o or t2,jv- The filled 
circles show H II regions with measured temperatures. The open circles show H II regions 
with measured t2,jv temperatures. The solid lines are the adopted cubic fits. 
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Fig. 3.— The [OIII]A5007)/H^ vs [NII]A6584)/H a diagram for subsamples A, B, and C and 
the total sample. T he solid line shows the d ividing line between Hll regions ionized by star 



clusters and AGNs (iKauffmann et al. 11200 



3). 
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Fig. 4. — The R3 - R2 diagram for subsumples A, B, and C and the total sample. The SDSS 
objects are shown by gray (light-blue in the electronic version) circles . The Hll regions in 
nearby galaxies from the compilation of iPilyugin. Vflchez. fc Contini I (120041 ) are shown by 
black triangles. 



-33 - 




Fig. 5. — The 0/H - N/0 diagram for subsumples A, B, and C. The SDSS objects are 
shown by gray (light-blue in the electronic versi on) open circles . The H II regions in ne arby 



galaxies with measured electron temperatures (IBresolin 1 120071 ; iBresolin et al. 
shown by black filled circles. 



2009b|) are 



-34- 



w 1.0 - 

CD 
-H 

X 0.9 - 

rd 
i — I 

rd 0.8 - 

tn 

U-1 0.7 - 

O 

U 0.6 - 

e o-s - 

3 

C 0.4 - 

TS 

g 0.3 - 

•H 

e 

S 0.1 - 



8.0 



n i i | i i i i | i i i i | i i i 

r-pa-- A ■ 



_l I I L 



j i i i_ 



10 

log(M s ) 



11 



T" 



T 



8.2 



8.4 8.6 

12+log(0/H) 



8.8 



12 



9.0 



■ i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i 




w 1.0 - 

Q) 
-H 

X 0.9 - 

rC 
i — I 

tfl 0.8 - 

tn 

M-i 0.7 |- 

O 

U 0.6 - 

CD 

3 

G 0.4 h 

S 0.3 h 

■iH 

0-2 h 

e 
a 



0.00 0.05 0.10 0.15 0.20 0.25 0.30 0.35 
redshif t 
i — i — i — i — | — i — i — i — i — | — i — i — i — i — | — i — i — i — i — | — i — i — i - 



0.0 




6.0 



7.0 7.5 

12+log (N/H) 



8.0 



8.5 



Fig. 6. — Normalized histograms of observed and derived properties for SDSS galaxies in 
subsamples A (short- dashed line), B (solid line), and C (long-dashed line). 
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Fig. 7. — The redshift - galaxy stellar mass diagram for subsamples A, B, and C. 
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Fig. 8. — The oxygen abundance - galaxy stellar mass diagrams for subsample B at 2=0.05 
(upper panel) and z=0.25 (lower panel). The Z-M5-O/H relation (Eqj9]) is shown by the 
solid line for z=0. 05 and by the dashed line for z=0.25. 
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Fig. 9. — The oxygen abundance - galaxy stellar mass diagrams for subsample C at z=0.05 
(upper panel) and z=0.25 (lower panel). The Z-Mg-O/H relation (Eq JTQj) is shown by the 
solid line for z=0. 05 and by the dashed line for z—0. 25. 
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Fig. 10. — Normalized histograms of oxygen abundance deviations from the Z-M^-O/H 
relation (upper panel) and nitrogen abundances deviations from the Z-M5-N/H relation 
(lower panel). In both panels, subsample B is shown by the solid line and subsample C by 
the dashed line. 
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Fig. 11. — The nitrogen abundance - galaxy stellar mass diagrams for subsample B at 
2=0.05 (upper panel) and z=0.25 (lower panel). The Z-M5-N/H relation (Eq JTTT) is shown 
by the solid line for z=0. 05 and by the dashed line for z=0.25. 
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Fig. 12. — The nitrogen abundance - galaxy stellar mass diagrams for subsample C at 
2=0.05 (upper panel) and z=0.25 (lower panel). The Z-M5-N/H relation (Eq JT2|) . is shown 
by the solid line for z=0. 05 and by the dashed line for z=0.25. 
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Fig. 13. — The oxygen and nitrogen abundances as a function of redshift for galaxies with 
masses 1O 1O -°M -1O 1O - 3 M (upper left panel) and 1O 1L2 M -1O 1L5 M (upper right panel). 
Oxygen abundances are shown by dark (black in the electronic version) filled circles, and 
nitrogen abundances by gray (light-blue in the electronic version) filled circles. The nitrogen- 
to-oxygen abundance ratios as a function of redshift for those galaxies are shown in the 
corresponding lower panels. 
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Fig. 14. — The oxygen abundance - galaxy stellar mass diagrams for the "all star-forming 
galaxies" (SFG) subsample at z=0.05 (upper panel) and z=0.25 (lower panel). The z-M s - 
O/H relation derived for subsample C (EqJTOl) is shown by the solid line for z=0. 05, and by 
the dashed line for 2=0.25. 
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Fig. 15. — The nitrogen abundance - galaxy stellar mass diagrams for the "all star-forming 
galaxies" (SFG) subsample at z=0.05 (upper panel) and z=0.25 (lower panel). The z-M s - 
N/H relation derived for subsample C (EqJT2l). is shown by the solid line for z=0. 05, and by 
the dashed line for £=0.25. 



